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Abstract We present two BeppoSAX observations of BL Lac (2200+420) as part of a multiwavelength campaign performed 
in 2000. The source was in different states of activity: in June, the X-ray spectrum was faint and hard (a ~ 0.8), with positive 
residuals towards low energies. In October we detected the highest [2-10] keV flux ever measured for the source. During this 
observation, the spectrum was soft (a ~ 1.56) up to 10 keV, while above this energy a hard component was dominating. The 
BeppoSAX data are confirmed by simultaneous RXTE short observations. During the first observation the soft X-ray flux was 
variable on timescales of a few hours, while the hard X-ray flux was almost constant. During the second observation, instead, 
the soft spectrum displayed an erratic behaviour with large variations (up to factors 3-4) on timescales smaller than 2 hrs. The 
analysis of the multiwavelength SED of October evidenced an intriguing feature: the optical and X-ray sections of the SED 
are misaligned, while in the prevailing standard picture, they are both thought to be produced via synchrotron emission. We 
suggested four scenarios to account for this discrepancy: a higher than galactic dust-to-gas ratio towards the source, the first 
detection of bulk Compton emission in the X-ray band, the presence of two synchrotron emitting regions located at different 
distances from the nucleus, the detection of a Klein-Nishina effect on the synchrotron spectrum. We evidenced the favorable 
and critical points of each scenario, but, at present, we cannot discriminate between them. 

Key words. BL Lacertae objects: general - X-rays: galaxies - BL Lacertae objects: individual: BL Lacertae (2200-1-420) 



1 . Introduction 

Blazars are radio-loud Active Galactic Nuclei producing vari- 
able non-thermal radiation in relativistic jets oriented close to 
the line of sight: the emission is therefore beamed and Doppler 
boosted (Blandford & Rees, 1978). They are characterized by a 
Spectral Energy Distribution (SED) displaying two broad fea- 
tures: the first, extending from radio to UV/X-ray is usually 
ascribed to synchrotron emission; the second, ranging from X- 
ray to y-ray, sometimes up to TeV energies, is attributed to 
inverse Compton scattering of seed photons by the same pop- 
ulation of synchrotron emitting electrons. The seed radiation 
field could be constituted by the synchrotron photons them- 
selves (SSC model, Maraschi, GhiselHni & Celotti, 1992) or 
by external photons, produced by the accretion disk (Dermer 
& Schlickeiser, 1993), by the Broad Line Region (Sikora, 
Begelmann & Rees, 1994; Ghisellini & Madau, 1996) or by hot 
dust (Blazejowski et al., 2000; Arbeiter et al., 2002). Different 
contributions of these fields can explain the observed blazar 
spectra. The BL Lac subclass is characterized by the absence 
or weakness of broad emission features and is divided in HBL 
and LBL (High and Low energy peaked BL Lacs) according to 
the radio to X-ray flux ratio (Padovani & Giommi, 1995). 
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BL Lac itself (lES 2200H-420) has been classified as a LBL on 
the basis of its radio-to-X broad band spectral index a„ = 0.85 
(Sambruna et al., 1996). It was first identified as the opti- 
cal counterpart of the radio source VRO 42.22.01 by Schmitt 
(1968); the presence of weak naiTow emission lines in its spec- 
trum allowed an accurate measurement of the redshift z - 
0.069 of the host elliptical galaxy (Miller & Hawley, 1977; 
Miller et al., 1978). It shows superluminal motions on m.a.s. 
scale(y3„pp ~ 3 -4 Mutel et aL, \99Q;l3app ~ 2.2-5.0/i"\ 
Denn et al., 2000). In spite of the definition of BL Lac objects 
as having featureless continua, in 1995, a survey performed 
by Vermeulen et al. (1995) revealed the presence of an Ha 
emission line with equivalent width of 7 , confirmed by sub- 
sequent observations (Corbett et al., 1996). During the same 
year EGRET, aboard the Compton Gamma Ray Observatory, 
detected a 4.4 cr excess above 100 MeV from its direction 
(Catanese et al., 1997). In the summer of 1997, BL Lac entered 
an exceptional flaring state with the highest X-ray flux ever 
recorded (Sambruna et al., 1999; Madejski et al., 1999) and a 
y-ray flux 4 times higher than in 1995 (Bloom et al., 1998) . 
In the X-ray band, BL Lac has been detected for the first time in 
1980 by the IPC (0.1-4 keV) and the MPC (2-10 keV) aboard 
the Einstein Observatory (Bregman et al., 1990). Since then 
the source has been observed many times by different satellites 
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such as EXOSAT (Bregman et al., 1990), GINGA (Kawai et 
al., 1991), ROSAT (Urry et al., 1996; Madejski et al., 1999), 
ASCA (Sambruna et al., 1999; Madejski et al., 1999), RXTE 
(Madejski et al., 1999) and finally BeppoSAX, which observed 
it in 1997 (Padovani et al., 2001) and twice in 1999 (Ravasio et 
al., 2002). 

During the second half of 2000, from July to December, the 
source has been the object of an intensive multiwavelength 
campaign (Bottcher et al., in prep.) which included two X-ray 
observations performed by BeppoSAX and was supplemented 
by a continuous long-term monitoring program by the Rossi 
X-ray Timing Explorer (RXTE), with 3 short pointings per 
week (Marscher et al., in prep.). 

During this campaign BL Lac has been observed in the radio 
band by the telescopes of the University of Michigan and of 
the Metsahovi Radio Observatory, while in the optical band it 
has been observed almost continuously by 24 telescopes in the 
context of an extensive WEBT campaign (Villata et al., 2002). 
Finally, HEGRA set an upper limit of 25% of the Crab flux 
above 0.7 TeV, after having accumulated a total of 10.5 h of 
on-source time during the autumn of 2000 (Mang et al., 2001). 
In this paper we will analyze in detail the BeppoSAX data of 
this campaign, comparing them with RXTE simultaneous ones 
and discussing the spectral and temporal behaviour of BL Lac 
in the X-ray band and in the whole radio-to-Te V energy range. 

2. BeppoSAX observations 

Thanks to its uniquely wide energy range (0.1-200 keV), the 
italian-dutch satellite BeppoSAX represents an ideal experi- 
ment for looking at blazars and expecially objects such as BL 
Lac, since it can detect the transition between the synchrotron 
and inverse Compton components of the SED (Tagliaferri et 
al., 2000; Ravasio et al., 2002). Therefore it allows to com- 
pare the simultaneous behaviour of extremely different parts of 
the emitting electron distribution. Boella et al. (1997 and refer- 
ences therein) report an extensive summary of the mission. 
BeppoSAX observed BL Lac (lES 2200H-420) twice during 
2000, since the July 26-27 measurements were soon inter- 
rupted (on-source time - 1.69 x 10"* s). Therefore we were 
given a second chance and a new observation started in October 
3 1 lasting until November 2, with a duration of 2.49 x 10"* s. In 
Table n] we report the exposures and the mean count rates for 
each BeppoSAX instrument. 

We performed our analysis on linearized and cleaned event files 
available at the BeppoSAX Science Data Center (SDC) on- 
line archive (Giommi & Fiore, 1998) using the software con- 
tained in the FTOOLS Package (XIMAGE 2.63c, XSELECT 
1.4b, XSPEC 10.00) and XRONOS 4.02. Data from MECS2 & 
MECS3 were merged by the SDC team in a single event file. 
Using XSELECT we extracted spectra and light curves from 
circular regions around the source of 8 and 4 arcmin radius for 
LECS and MECS, respectively. We extracted event files also 
from off-source circular regions, in order to monitor the back- 
ground behaviour during our measurements. Since the LECS 
and MECS backgrounds are not uniformly distributed across 
the detectors, after having checked the constancy of the ex- 
tracted background Ught curves, we choose to use the back- 



ground files obtained from long blank field exposures, avail- 
able at the SDC public ftp site (Fiore et al., 1999; Parmar et al., 
1999). 

2.1. Spectral analysis 

The spectral analysis was performed with XSPEC 10.0, using 
the updated response matrices and blank-sky background 
files (01/2000) available at the SDC public ftp site. The 
LECS/MECS normalization factor was fixed at 0.72, while 
the PDS/MECS was fixed at 0.9 to be consistent with our 
previous works (Ravasio et al., 2002): these values are within 
the acceptable range indicated by SDC (LECS: 0.67-1, PDS: 
0.77-0.93, Fiore et al., 1999). 



2.1.1. July 26-27 

During this observation, the source was not detected by PDS 
because of the short on-source time and because of the intrinsic 
weakness of the source itself. For the same reasons and because 
of the high galactic absorption, the detection was uncertain also 
at LECS low energies: therefore we proceeded to the analysis 
only in the [0.6-10] keV range, fitting the extracted spectrum 
with a single power law model. 

We repeated the procedure three times, letting the interstellar 
absorption parameter either free to vary, or fixed at the value 
Nh = 3.6 X 10^' cm"^. The latter is obtained adding the galactic 
value from 21 cm measurements Nh - (2.0 + 1.0) x 10^' cm~^ 
(Elvis et al., 1989) and the absorption due to the molecular 
clouds observed along the line of sight (Lucas & Liszt, 1993). 
Finally, we used the fixed value Nh - 2.5 x 10^' cm"^, to be 
consistent with previous works (Ravasio et al., 2002; Sambruna 
et al., 1999; Madejski et al., 1999). 

Letting the absorption parameter Nh free, the 0.6-10 keV spec- 
trum is well fitted by a single power law with spectral in- 
dex a = 0.75 + 0.15 ixlld.o.f. = 0.96/50). We obtained 
Nh = 2.0 + 1 x 10^' cm"^, the large uncertainties are due to 
the low count rate in the soft X-ray band. 
A single power law model with Nh fixed at 3.6 x 10^' cm"^ 
leaves large positive residuals below ~ 1 keV. This suggests 
that we are using an exceedingly high absorption column or 
that near the low boundary of our energy range we are detecting 
the transition between a soft and a hard spectral component. 
These residuals become less relevant adopting the intermedi- 
ate value Nh = 2.5 x 10^^ cm"^ (fig. 0: in this case a sin- 
gle power law model with ax - 0.8 + 0.1 fits the data well 
0^2/^.0./. = 0.96/51). 

According to an F-test, a broken power law model does not im- 
prove significantly the quality of the fit. The best fit parameters 
for each model are shown in table|2] which reports also the flux 
at 1 keV and the integrated flux in the 2-10 keV energy range. 

2.1 .2. October 30- November 2 

The October observation is more interesting because of the 
longer duration and the higher state of the source. This allows 
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LECS 




MECS 




PDS 


Date 


exposure 


count rate" 


exposure 


count rate* 


exposure 


count rate'^ 




(s) 


xlO"^ cts s"' 


(s) 


xlO"- cts s"' 


(s) 


xlO"^ cts s"' 


26-27 July 2000 


16917 


4.99 ± 0.20 


23309 


8.23 ±0.21 


1.04 X lO'' 


-1.20 ±5.73 


31 October- 














2 November 2000 


24899 


28.5 + 0.38 


33661 


33.5 ±0.33 


1.88 X 10^ 


12.92 ± 4.27 



Table 1. Journal of Be;?/?o5AX observations. 0.1-lOkeV; 1.5-lOkeV; 12-100 keV. 
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Xlld.o.f. 
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keV 




^Jy 
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cm"^ 


S-' 
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2.0 ± 1 


0.75 ±0.15 






1.1 


5.8 


X 10- 


12 


0.96/50 


fixed 


















3.6 


0.9 ±0.1 






1.3 


5.8 


X 10- 


12 


1.07/51 


fixed 


















2.5 


0.8 ±0.1 






1.1 


5.8 


X 10- 


12 


0.96/51 



31 October - 2 November 2000 



Nh 

xlO^' cm-2 


ai 


Eft Q'2 Fu-fV 

keV yuJy 


F2-10teV 

erg cm"- s"' 


xild.o.f. 


Power law models 


free 

3.0 ± 0.3 


1.63 ±0.06 


13.0 


2.1 X 10-" 


0.69/76 


fixed 
3.6 


1.71 ±0.03 


14.6 


2.1 X 10-" 


0.96/77 


fixed 
2.5 


1.56 ±0.03 


11.9 


2.1 X 10-" 


0.87/77 


Broken power law models 


fixed 
2.5 


1 45+"' 


2.2+^11 1.65 ±0.08 11.4 


2.1 X 10-" 


0.68/75 



Table 2. LECS + MECS spectrum best-fit parameters. 




13 5 10 

channel energy (keV) 



Figure 1. LECSh-MECS 26-27 July spectrum of BL Lac. The 
spectrum is modeled with a power law and a fixed Nh = 2.5 x 
10^' cm"^. We choose this value to be consistent with previous 
work (Ravasio et al., 2002; Madejski et al., 1999; Sambruna et 
al., 1999) 



us to analyze a wider spectral range, from 0.3 keV up to 50 
keV, thanks to the detection by the PDS. 
As before, we performed the fits letting Nh free, then fixing it to 
its maximum value and finally to Nh - 2.5 x 10^' cm"^. In the 
first case, the LECS + MECS spectrum was well fitted by a soft 
single power law model, with energy index ax - 1.63 + 0.06 
(Nh = 3.0 + 0.3 X 10^' cm"^). A similar result is obtained also 
keeping Nh fixed to 2.5 x 10^' cm"^, while the model with the 
highest absorption value leaves positive residuals increasing to- 
wards low energies. 

We repeated the procedure using a broken power law model 
and fixing Nh = 2.5 x 10^' cm"^, for consistence with previous 
works (e.g. Ravasio et al., 2002). The LECSh-MECS spectrum 
in the latter case is well fitted by a convex curve steepening be- 
yond ~ 2 keV. The best-fit parameters of each LECS + MECS 
spectral model are listed in Table|2] An F-test suggests that the 
addition of two parameters gives a 99.9% probability of im- 
proving the quaUty of the fit. 

Subsequently we included also PDS data in the analysis. A sin- 
gle power law model leaves large positive residuals towards 
high energies (fig.|2j: the [0.1-50] keV spectrum seems to be 
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Figure 2. LECS+MECS+PDS 31 October-2 November spec- 
trum of BL Lac. Both a single and a broken power law model 
cannot reproduce the whole observed spectrum. PDS data lie 
above the power-law best fit of the of LECS+MECS spectrum. 
The LECS+MECS spectrum is well fitted by a convex broken 
power law. 

concave. PDS data lie largely above the extrapolation of the 
LECS and MECS spectrum: this can be explained as the ris- 
ing of a second, hard spectral component. However, the en^or 
bars are too large and we cannot constrain a second compo- 
nent in the model. Thus, in order to have an idea of this hard 
component spectral shape, we fitted the MECS + PDS spectra 
with a broken power law model, keeping the low energy spec- 
tral index fixed to the value obtained from the fit of the LECS 
+ MECS spectrum (~ 1 .65, independently from the absorption 
parameter chosen). The best fit of the two extra parameters are 
Eb ~ 9.9 keV and q-pds ~ 0.56, but we are not able to estimate 
the uncertainties. 



els of RXTE spectra are consistent with those of the MECS, 
confirming the results obtained from BeppoSAX spectral anal- 
ysis. 

2.2. Temporal analysis 

We performed this analysis with the XRONOS 5.16 package, 
grouping data in different time bins to reduce uncertainties. 
During the observation of July 2000 the source was in a faint 
state of activity so we choose bins of 3600 s. During the sec- 
ond observation, the source emission was stronger so we could 
rebin our data in smaller bins of 600 s. In our analysis we ex- 
cluded the bins with less than 20% of effective exposure time, 
to reduce weighting errors. 

2.2.1. 26-27 July 

In order to have information on the behaviour of the source at 
different X-ray wavelenghts we extracted a soft ([0.7-2] keV, 
LECS) and a hard ([2-10] keV, MECS) X-ray light curve. The 
source displays large variability in the low energy band: the 
0.7-2 keV flux increases by a factor > 2 on time scales of ~ 4 
hours, then fades to previous values in ~ 5 - 6 hours (Fig.|3}- 
In the harder energy band, variability events are less evident 
but the source is still active, with flux variations of 25 - 30% 
(X^-test constancy probability ~ 10%). The two light curves 
are only marginally correlated: the null correlation probability 
is quite high (~ 25%). 

This behaviour was noticed also in the June 1999 BeppoSAX 
observation of BL Lac (Ravasio et al., 2002): in the high energy 
MECS band (4-10 keV) the source was not variable, while the 
[0.3-2] keV LECS and the [2-4] keV MECS light curves were 
displaying a flare of a factor ~ 4. 



2.1.3. Simultaneous RXTE observations 

In addition to BeppoSAX, the 2000 campaign was covered in 
the X-ray band also by Rossi X-ray Timing Explorer (Bradt, 
Rotschild & Swank, 1993), which provided three short expo- 
sures a week, covering the whole duration (Marscher et al, in 
prep.). We analyzed the two RXTE pointings temporally clos- 
est to the BeppoSAX ones. In July the two observations were 
exactly simultaneous while in November RXTE was lagging 
BeppoSAX only by 1 hour. Therefore we are given the oppor- 
tunity to test our results with a totally independent set of data. 
We restrict our analysis to PCA data (Jahoda et al., 1996), an in- 
strument composed by 5 passively collimated independent X- 
ray detectors sensitive to the 2-60 keV range. 
We choose to compare the [3-15] keV RXTE spectrum to the 
[1-10] keV MECS using a power law model with fixed absorp- 
tion parameter (Nh = 2.5 x 10^' cm"^). The model parameters 
are reported in Table |3l together with the log of the observa- 
tions. 

Both RXTE spectra are well fitted by power law models: dur- 
ing the observation of July, RXTE detected a hard spectrum 
(a - 0.9), while in October, it detected a soft component 
(a = 1 .45). The slope and the normahzation of the best-fit mod- 



2.2.2. 31 October-2 November 

During the October BeppoSAX observation, BL Lac displayed 
very fast and remarkable flux variations in the whole energy 
range covered by the LECS and the MECS (Fig.©. The light 
curves are characterized by a similar, erratic behaviour, with 
flux variations of factors ~ 3 - 4, even on timescales of ~ 1 
hour The three light curves are highly con^elated (P> 99.9%) 
even with a short temporal binning of 10 min. 
We analyzed also the PDS light curve, which turned out to be 
constant: a test gives a ~ 96% constancy probability. Because 
of the large errorbars, however, we would not be able to detect 
variations smaller than a factor of 3. For BL Lac, this X-ray be- 
haviour is not unprecedented: as already mentioned, during the 
June 1999 BeppoSAX observation, the [0.3-2] keV flux dou- 
bled in 20 min and faded to previous values in a similar time. In 
less than 2 hours a complete flare with a flux variation of a fac- 
tor ~ 4 was observed. The amplitude of this event was highly 
frequency-dependent: while the flare was extremely prominent 
in the [0.3-2] keV and in the [2^] keV curves, at higher fre- 
quencies the flux remained constant. 

The best way to characterize the temporal behaviour of the 
source would be the calculation of the power density spectrum: 
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26-27 July 


Instrument 


Time start 


Time end 


Duration 

(sec) 


a 


0"Jy) 


p. 

(xl0-'2) 




PCA 


26/7/00 


26/7/00 


2208 




1.4 


6.3 


0.42/25 


3-15 keV 


18:23:44 


19:00:32 










MECS 


26/7/2000 


27/7/2000 


23309 


0.8 + 0.1 


1.18 


5.8 


0.86/43 


1-10 keV 


10:12:39 


06:43:33 












30 October-2 November 


PCA 


2/11/2000 


2/11/2000 


2000 


1 45+"" 


10.3 


19.7 


0.45/25 


3-15 keV 


10:56:16 


11:29:36 










MECS 


31/10/2000 


2/11/2000 


33661 


1.6 ±0.05 


12.7 


19.7 


0.61/58 


1-10 keV 


20:46:55 


09:59:28 













Table 3. ': erg cm ^ s RXTE and BeppoSAX simultaneous observations log and spectral fit parameters. We adopted a power 
law model with fixed absorption parameter: Nh = 2.5 x 10^' cm"^. 



26-S7/7/S000 



Bin lime = 3600 s 




15 20 25 

Time (hours) 



31/10-2/11/2000 



Bin time=600 s 



0.4 r 
0.3 r \ 



I 

£0.2 -X I * i 



+ 



0.3 



1 0.2 i 



1 ij^ S^lii 



0.15 



Pill 



-I 1 1 1 1 



30 40 

Time(hours) 



Figures. LEGS 0.7-2 keV (top panel), MECS 2-10 keV (bot- 
tom panel) light curves of BL Lac during the July 2000 
BeppoSAX observation. Note that in the LECS range the 
source displays larger variability than at higher energies. 



Figure 4. LECS 0.5-2 keV (top panel), MECS 2-4 keV (mid 
panel) and MECS 4-10 keV (bottom panel) light curves of BL 
Lac during the October-November 2000 BeppoSAX observa- 
tion. 



however this tool is not appropriate for unevenly sampled light 
curves like BeppoSAX ones. The large observational gaps, the 
short duration of the run and the limited statistics caused by the 
faintness of the source also makes the use of alternative tech- 
niques very difficult, such as the discrete correlation function 
(Edelson & Krolik, 1988) or the structure function calculation 
(Simonetti, Cordes & Heeschen, 1985). 

We can still characterize the variability degree of the source 
using two common estimators: the normalized excess variance 
parameter cr^^^ and the minimum doubling timescale short 
(Zhang et al., 1999; Fossati et al., 2000). The first parameter 
is defined as the normalized difference between the variance of 
the light curve and the variance due to measurement errors: it 
quantifies the mean variability of the source. The second pa- 
rameter, instead, represent a measure of the fastest significant 



timescale of the source (Edelson, 1992). Assuming that each 
point of the light curve is described as (f,,./i) we define the "dou- 
bling time" as 



Tij - 



(1) 



where Ar,y = f,- - tj, Afij = - fj and fij = (fi + fj)/2. 
For each of the discussed energy ranges ([0.5-2] keV-LECS; 
[2^] keV-MECS; [4-10] keV-MECS) we produced 4 differ- 
ent light curves with time bin of 500, 1000, 1500 and 2000 sec, 
respectively. For each curve we computed the cr^^^^ and the av- 
erage of the 5 smallest Tjj (calculated for each data pair (i,j) of 
the light curve) with fractional error lower than 25%. We define 
Tshort as their weighted average. In Table|3]we summarize our 
results. 
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30/10-2/11 2000 



Energy band 


<Count rate>° 




< o-'Ls >" 




(keV) 


(10-2 


(10-2) 


(10-2) 


(ksT 


LEGS 


17.6 


15.2 


12.8 


6.3 


0.5-2 keV 


±7.2 


±3.3 


±1.6 


±0.6 


MECS 


14.7 


11.4 


10.1 


5.5 


2-4 keV 


±5.3 


±2.5 


±1.5 


±0.5 


MECS 


7.4 


7.9 


7.5 


7.1 


4-10 keV 


±2.5 


±2.1 


±1.1 


±1.3 



Table 4. Variability parameters. ") From 600 sec. rebinned light curves. Weighted mean of data taken from 500, 1000, 1500, 
2000 sec. rebinned light curves. 



The highest frequency light curve is the least variable, as can 
be noticed also by looking at Fig. |3 probably this can be ex- 
plained partially as an effect of the larger measurement error 
caused by poorer statistics and partially as a lower variability 
at higher energies. 

Within the relatively large errors, no significant differences 
are found between the Tshori values of the three light curves: 
the source seems to have a characteristic minimum variability 
timescale of ~ 1.5-2 hrs. 

3. Spectral Energy Distributions 

After having performed the spectral analysis we are now able 
to build the SEDs of BL Lac relative to the BeppoSAX 2000 
observations and to compare them with other historical multi- 
wavelength SEDs. 

As we have discussed above, during 2000 BeppoSAiX detected 
the source in two completely different states of activity: in July 
the source was detected only in the [0.6-10] keV range, dis- 
playing a faint (F2-10 = 5.8 x 10"'^ erg cm-^ s-') hard spec- 
trum. In October-November, instead, BL Lac was displaying a 
very intense (F2-10 = 2.1 x lO"" erg cm-^ s-') soft spectrum 
up to ~ 10 keV, while at higher energies a hard component, de- 
tected by the PDS up to 45 keV, was dominant. 
Also in the optical band the source was in different states: the 
increase of the optical flux, in fact, was the reason for prolong- 
ing the multiwavelength campaign beyond August 2000. The 
optical fluxes measured simultaneously to the BeppoSAX ob- 
servations differed by 40% between the two pointings: in the 
core of the campaign (26-27 July) the source average R-band 
magnitude was = 14.08, while during the second X-ray run 
it was m^ - 13.74 (Viflata et al., 2002). After having dered- 
dened the data using Ab = 1.42 (Schlegel et al., 1998), we 
calculated the optical spectral indices using weight-averaged 
UBVRI fluxes (see values in fig.|5lcaption): when fainter, the 
spectrum is softer (ffop, ~ 1.84+0.01), while during the autumn 
it is harder (oopt ~ 1.58 + 0.04). 

In Fig. Is] we plot the two simultaneous multiwavelength SEDs 
of BL Lac. We report also the upper limit by HEGRA, 7.7 x 
10-'^ photons cm-2 s-' above 0.7 TeV (Mang et al., 2001), 
which is 25% of the Crab level: this means Fo vxev - 5.13 x 
10"'^ X (ffTev) Jy, where ffxev is the energy index of the TeV 
spectrum (assumed to be = 1.5 in the plot). 



3.1. Historical SEDs 

As mentioned in the introduction, blazar SEDs are character- 
ized by two broad humps extending, the first, from radio to 
UV/X-rays and the second from X-rays to y-ray or even TeV 
energies. The first component is usually interpreted as syn- 
chrotron emission from high-energy electrons in relativistic 
motion, while the second is explained as inverse Compton scat- 
tering of seed photons of still debated origin. 
The multiwavelength history of BL Lac follows that of its X- 
ray observations which can be traced back to 1980, the year 
of the first X-ray detection, carried out by Einstein (Bregman 
et al., 1990). BL Lac was then observed also by EXOSAT 
(Bregman et al., 1990), GINGA (Kawai et al., 1991), ROSAT 
(Urry et al., 1996; Madejski et al., 1999), ASCA (Sambruna et 
al., 1999; Madejski et al., 1999), RXTE (Madejski et al., 1999) 
and by BeppoSAX (Padovani et al., 2001 ; Ravasio et al., 2002). 
In table |5l we report the published spectral parameters for each 
X-ray observation, together with simultaneous optical ones. 
The X-ray spectra of BL Lac exhibit a variety of different 
shapes: sometimes they are concave, hardening towards high 
energies, and sometimes they are hard in the whole X-ray 
band. In 1980, IPC ([0.1-4] keV) and MPC ([2-10] keV) aboard 
Einstein detected quite different spectral indices (ai = 2.2; 
a2 - 0.8; Bregman et al., 1990) and a similar shape was de- 
tected also by ASCA in November 1995 (Sambruna et al., 
1999) and by BeppoSAX in June 1999 (Ravasio et al., 2002), 
as well as during this campaign. In other epochs the X-ray 
spectrum of BL Lac was displaying only a hard component. 
In June and July 1988 GINGA detected faint X-ray spectra 
which can be attributed to inverse Compton emission in a lep- 
tonic jet model (Kawai et al., 1991), similarly to the obser- 
vation of ROSAT (1982), to those performed by BeppoSAX 
in November 1997 (Padovani et al., 2001), December 1999 
(Ravasio et al., 2002) and to our observation of July 2000. All 
the hard X-ray spectra are similar to each other except for one 
case: during July 1997 BL Lac was in a high state, consistent 
with Compton emission, with a 2-10 keV flux ~ 3 - 4 times 
higher than that observed in 1995 by ASCA or in June 1999 
by BeppoSAX; the high state was clearly detected also in the 
y-ray by EGRET (Bloom et al., 1997). The complexity of the 
behaviour of the source, suggested by the data reported in table 
|5] is further highlighted by the comparison of the multiwave- 
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Figures. Left panel: 26-27 July SED. Right panel: 31 October-2 November SED. The black butterflies represent RXTE si- 
multaneous data. The HEGRA upper limit is calculated assuming a spectral energy index axev = 1-5 above 0.7 TeV. The op- 
tical slopes are calculated on the weigthed average UBVRI fluxes (Ab = 1.42). July fluxes (mJy): < Fi >= 5.17 + 0.07; < 
Fr >= 8.12 ± 0.01; < Fv >= 12.08 + 0.04; < Fb >= 16.54 + 0.02; < Fu >= 22.28 + 0.06. October fluxes (mJy): 
<Fi >= 8.26 ± 0.09; < Fr 11.96 + 0.05; < Fy >= 17.56 + 0.06; < Fb 22.43 + 0.08; < Fu >= 27.97 ±0.18. 



Date 




0"Jy) 


F2-mkeV 

(10-'2 c.g.s.) 




Fv 
(mJy) 


June 1980 


2.2 ± 1.0 


1.6 


7.1* 


2.2 ± 0.2 


5.75 


December 1983 


(2.2) 


0.8 




2.1 ±0.4 


9.55 


June 1988 


0.71 ±0.07 


1.27 


7.7* 


2.2 ±0.1 


4.2 


July 1988 


1.16 ±0.24 


1.35 


4.1* 


3.21 ±0.02 


6.46 


December 1992 


0.94 ± 0.46 


0.88 


3.7* 






November 1995 


0.96 ± 0.04 


2.22 


8.9 


0.8 ± 0.2 


5.7 


July 1997 


0.59 ± 0.03 


3.55 


26 


1.85 ±0.4 


27.9 


November 1997 


0.89 ±0.13 


2.4 


11.2 






June 1999 


1.57 ±0.25 


0.77 


6.5 


1.2 ±0.2 


26.4 


December 1999 


0.63 ± 0.06 


0.96 


12.3 


1.8 ±0.3 


18.6 


July 2000 


0.8 ±0.1 
0.8 ±0.1 


1.1 
1.1 


5.8 
5.8 


1.84 ±0.01" 
0.94 ± 0.01' 


12.08" 
25.60* 


October 
November 2000 


1 45+" ' 
1 45+" ' 


11.4 
11.4 


21 
21 


1.58 ±0.04" 
0.64 ± 0.02* 


17.56" 
37.23* 



Table 5. June 1980 Einstein data: Bregman et al., 1980. December 1983 EXOSAT data: Bregman et al., 1990. June 1988-July 
1988 GINGA data: Kawai et al., 1991. December 1992 ROSAT data: UiTy et al., 1992. November 1995 ASCA data: Sambruna 
et al., 1999. July 1997 RXTE data: Madejski et al., 1999. November 1997 BeppoSAX data: Padovani et al., 2001. June and 
December 1999 BeppoSAX data: Ravasio et al., 2002. *: from our calculations. using Ab - 1.42. *: using Ab = 2.5: see 
discussion. 



length SEDs: in fig.|6lwe plot the results of all the multiwave- 
length campaigns performed on BL Lac. 
As evidenced in fig.|5l&|6land in table |5] while during July 
2000 we observed the source in a normal state, in the end 
of October BL Lac was very active. During this run, in fact, 
BeppoSAX detected the highest soft X-ray flux and an inte- 
grated [2-10] keV flux which is only slightly smaller than that 



of the flare of July 1997, when BL Lac was displaying a hard 
X-ray spectrum. 



4. Discussion 

X-ray observations are crucial for blazar study because at these 
energies we can often observe the transition between the two 
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Figure 6. Simultaneous multiwavelength SEDs of BL Lac. Left-top panel: July and October-November 2000 SEDs 
(BeppoSAX). We have reported optical data corrected for Ab = 1.42, to be consistent with the other historical observations. 
Right-top panel: November 1995 (ASCA) and July 1997 (RXTE) SEDs. Left-bottom panel: June and December 1999 SEDs 
(BeppoSAX). Right-bottom panel: June 1980 (Einstein), December 1983 (EXOSAT), June and July 1988 (GINGA) multiwave- 
length SEDs. The dotted butterfly represent the December 1982 ROSAT data. 



emission mechanisms. The two X-ray observations that we 
presented here clearly show this: 

- 26-27 July: during this observation BeppoSAX detected a 
hard spectrum in the LECS-MECS range, with large pos- 
itive residuals below 1 keV. The hard component can be 
attributed to inverse Compton emission while the residu- 
als can be explained as the very soft tail of the synchrotron 
component (the data however does not allow us to set lim- 
its on the spectral index of this component neither to asses 
its consistency with the optical data). This picture is sup- 
ported by the temporal behaviour of the source: the [2-10] 
keV light curve represents the inverse Compton spectrum 
produced by the less energetic particles and is therefore al- 
most constant on the timescale of a single observation. The 
[0.7-2] keV curve, instead, is influenced also by the soft 
synchrotron emission produced by highly energetic elec- 
trons which can account for the observed variability. 



- 31 October - 2 November: also in this observation 
BeppoSAX detected the transition between two emission 
mechanisms. The LECSh-MECS spectrum is soft, while 
the PDS observed a hard component. In the usual scenario, 
the hard component is attributed to inverse Compton emis- 
sion, while the soft spectrum is explained as synchrotron 
emission by very energetic electrons (y ~ 2.5 x 10"^ - 10^ 
(B/1G)"'''^(5/10)"'^^) which have a very short coohng 
time (t.^,„ -3x10^- lO\B/lG)-^'\6/W^ s). This could 
also explain the very fast variability detected in the whole 
LECS-MECS energy range. 

This scenario, however, could be inadequate. Looking 
at fig. 13 one can notice immediately a strange feature 
in the SED of October-November: the X-ray data lie 
above the extrapolation of the optical spectrum, while 
both of them should be produced by the same emission 
mechanism, the synchrotron. Even correcting for the host 
galaxy contribution (which is almost negligible for the 
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level of activity of BL Lac during our observations) we are 
not able to reconcile X-ray and optical data. However, we 
can conceive at least four scenarios which under certain 
conditions can explain the observed spectral "glitch". They 
are: I) a variable local absorbing column along the line 
of sight; II) bulk Compton radiation; III) two different 
synchrotron emitting regions; IV) Klein-Nishina effect on 
the synchrotron spectrum . 

CM 

E 

o 
w 

4.1. Differential absorption ^ 

From the BeppoSAX X-ray analysis we had Nh = 2.5 x 10^'^ 
cm"^, consistent with previous observations; using the dust-S* 
to-gas ratio suggested by Ryter (1996) 

Av = 4.5 X 10"^^ Nh cm"^ (2) 

and using E(B-V)=0.329 as indicated by the NED, the corre- 
sponding optical absorption is Ab = 1.45, very similar to the 
value reported by Schlegel et al. (1998), Ab = 1.42. However, 
using these values, the X-ray spectrum lies largely above the 
optical spectrum extrapolation, as illustrated in Fig. Even 
correcting the optical data for the host galaxy contribution the 
two sections do not connect continuously. In Fig.0we plotted 
the R and B band BL Lac fluxes corrected using R;,oj, = 15.55 
and B/,oj/ = 17.15 (Villata et al., 2002): the misalignement is 
even worsened. 

If we try to account for this discrepancy using different absorp- 
tion values, e.g. Ab - 2.5, and consequently Nh = 4.8 x 10^' 
cm"-, we still obtain two spectral sections which can not be 
smoothly connected (Fig.Q. 

Therefore we fixed the X-ray absorption to Nh = 2.5 x 10"' 
cm"^ and varied the optical absorption. 

Using Ab ~ 4.2, U-filter data lie on the X-ray extrapolation, 
but the optical spectrum would be very hard and the connec- 
tion discontinuous. Using an intermediate value, Ab = 2.5, the 
optical and X-ray spectra connect continuously on a parabolic 
curve, evidenced in Fig.0 With this absorption value, we ob- 
tain a hard optical spectrum (ofy,,/ = 0.94 + 0.01 and aoct = 
0.64 + 0.02, for the summer and the autumn observations re- 
spectively) and consequently higher synchrotron peak frequen- 
cies. A dust-to-gas ratio higher than that reported by Ryter 
(1996) could therefore account for the observed optical to X- 
ray misalignement. 

We can compare this behaviour with other multiwavelength 
campaigns: in Fig. |6l we plot all the published SEDs of BL 
Lacertae. Only in the SED of 1980 and in that of June 1999 
we observe synchrotron emission both in the optical and in 
the X-ray band. In all the other cases, only the Compton 
component was detected, thus one can not see any misaligne- 
ment between the optical and X-ray synchrotron components. 
While the BeppoSAX synchrotron spectrum of June 1999 con- 
nect continuously with the simultaneous optical data (using the 
galactic dust-to-gas ratio; Ravasio et al., 2002), in the 1980 
SED the optic al-UV data seem to lie below the extrapolation 
of the Einstein IPC spectrum (Bregman et al., 1990), recreating 
the misalignement observed in October-November 2000. 
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Figure?. Optical to X-ray SED of October-November 2000 
BeppoSAX observation. Filled circles represent data corrected 
for Ab = 1.42 corresponding to Nh - 2.5 x 10^' cm"^. Filled 
triangles represent data corrected for Ab = 2.5 = Nh = 4.8 x 
10"' cm"^. Open circles represent raw optical data and optical 
data dereddened with Ab = 4.2. The two stars represent the B 
and R-band data dereddened with Ab = 1 .42 and corrected for 
the host galaxy contribution. The two dotted lines represent the 
extrapolation of the BeppoSAX spectrum to the optical band 
and a parabolic curve which connects X-ray and optical data 
dereddened with Ab - 2.5. See text for discussion. 



During the 1980 IPC observation, in fact, BL Lac was dis- 
playing a very steep synchrotron spectrum in the [0.1^] keV 
range (a - 2.2 + 1.0, using however a very high absorption 
parameter Nh - 1.25 x 10^^ cm"^; Bregman et al., 1990), 
hardly connectable with the very soft simultaneous optical- 
UV spectrum. The Einstein data, however, are affected by very 
large uncertainties, which make difficult to determine the exact 
shape of the SED: this is evidenced by the results published by 
Worrall & Wilkes (1990), which reported an IPC spectral in- 
dex a - 1.34^^ 3. The uncertainties on the IPC spectral index 
are such that no firm conclusion about the reality of the optical- 
to-X-ray misalignement in the 1980 SED can be drawn. 
The strange optical/X-ray misalignement observed in 2000 
(and maybe in 1980) is not detected in the only other 
multiwavelength campaign that shows synchrotron emission 
in the optical and X-ray bands (June 1999, BeppoSAX). 
Nevertheless, we can rely on the goodness of our data: the 
BeppoSAX spectrum is confirmed by simultaneous RXTE 
data, as described in the previous sections, while the optical 
data are confirmed by different observatories (Villata et al., 
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2002). The optical/X-ray misalignement is therefore real. The 
absorption along the line of sight can not account for it, unless 
we assume a sudden, large increase of the dust-to-gas ratio in 
the interstellar material. 

This hypothesis, however, is not unlikely, since the line of sight 
towards BL Lac is partially covered by a low surface bright- 
ness interstellar nebulosity which is very variable (Penston & 
Penston, 1973; Sillanpaa, priv. comm.). If these clouds are 
dusty, they could account for the dust-to-gas ratio excess 
needed to reconcile the optical and X-ray data in the SED of 
October-November 2000 (and possibly also in the 1980 SED; 
Sillanpaa et al., 1993). Their proper motion could indeed ex- 
plain the misalignement in the SED of autumn 2000 and its 
absence in that of June 1999. 

4.2. Bulk Compton radiation 

Alternatively, the observed misalignement could be claimed as 
the first detection of the so called "bulk Compton" emission. 
This possibility was predicted by Sikora et al. (1997): he postu- 
lated the existence of a population of "cold" electrons in addi- 
tion to the relativistic ones producing synchrotron and inverse 
Compton emission. These cold electrons, however, have bulk 
relativistic motion with respect to the radiation fields produced 
by the accretion disk and by the broad line region and can in- 
verse Compton scatter these photons up to frequencies 



vbc 



(3) 



where vq is a characteristic frequency of the external radiation 
field and F is the bulk Lorentz factor of the jet. If the seed ra- 
diation peaks in the UV band, as the accretion disk and the 
broad Une region emission do, the bulk Compton component 
will peak in the soft X-ray band accounting for the misaligne- 
ment we observe in the 2000 BL Lac SED. 
Assuming a conical jet with opening angle Oj ~ l/T and the 
conservation of the flux of electrons along the jet (n^ cc l/r^), 
Sikora et al.(1997) estimated the observed amount of bulk 
Compton radiation Lsc- 



--BC 



2r^ 



ne(rmm)rminO-T 



UV 



(4) 



where r^in is the distance from the apex of the cone at which 
the jet is fully developed, Luv is the external radiation field and 
( - ^(r) is the fraction of the UV radiation contributing to the 
isotropic field at a distance r from the jet apex (Sikora et al., 
1997). 

In October-November the X-ray luminosity was Lx - 10*^ 
erg s"^. Assuming that the bulk Compton emission peaks at 
V = 10'^ Hz and vq = 10'^ Hz, from Eq.|3]we have ~ 10. 
Since Vermeulen et al. (1995) observed broad Ha and H/3 emis- 
sion lines (confirmed by Corbett et al., 2000), from their data 
we can evaluate Luv in an indirect way. Postulating a fixed line 
ratio and following the method described in Celotti, Padovani 
& Ghisellini (1997) we obtain an average value Lblr ~ 5 x 10"*^ 
erg s If the Broad Line Region covering factor is ~ 10%, the 



Using these values and assuming Lbc = Lx we obtain the par- 
ticle number density in the observer frame 



^f(^min) 



6x 10^ 



Tjo^O.! '"min,15 



cm 



-3 



(5) 



where Fio = F/10, ^o.i = f/0.1 and r^i^js = r^in/lO'^cm. 
This number particle density, necessary to produce the ob- 
served X-ray spectrum, puts constraints on the jet composition. 
In fact, if we suppose a one-to-one proton-electron plasma, we 
can easily evaluate the total kinetic power of the jet (Celotti & 
Fabian, 1993): 



^KIN ~ ^(fmmdj) n'^MpC V C ~ 
, „47 ^lo.l _i 

~ 10 — erg s ' 

Fio (o.i 



(6) 



where - «p/F is the proton number density in the jet comov- 
ing frame. This value exceeds the Lkin estimated by Celotti, 
Padovani & Ghisellini (1997) by a factor ~ 30. To produce the 
observed X-ray spectrum via bulk Compton mechanism and to 
be consistent with the result of Celotti, Padovani & GhiselUni 
(1997) the proton-electron ratio must be smaller: the jet must 
be pair loaded. 

This is not forbidden since the optical thickness for Thomson 
scattering (T(r) = neO-rOjr) at rmin is 



T(rmm) = 



2 



(7) 



0.1 



Since from our calculations T(rniin) < 1 (and decreasing further 
out) and the expansion timescale at r^in (fg^p ~ 8 x lO'^rn^n.isrio 
s) is smaller than the annihilation timescale (?«„„ ~ 2 x 
10^ F|Q ^o.i'"min,i5 s) the pairs will survive along the jet. 
After having ruled out the hypothesis of a proton-electron jet 
and having requested the presence of pairs in the jet, it is now 
fundamental to understand if the particle number necessary to 
produce the bulk Compton emission is sufficient to produce the 
observed SEDs. The particle number conservation we have just 
demonstrated implies that «e oc l/n. Therefore, at the distance 
where the radiation is emitted, tsed ~ 10'* - 10'^ cm, the par- 
ticle number density would be 



"('"SEd) = n('"min)(-^^) ~ 
V '"SED / 

6 X 10 ^ ^ cm 

Tjo ^0.1 ''SED.IV 



(8) 



disk luminosity is Li 



UV 



5 X lO"*"* erg s 



This value is greater than what we needed to model old BL Lac 
seds: ~ lO^/ymin cm""^. Only a fraction of the electrons car- 
ried by the jet need to be accelerated in order to produce the 
observed radio to y-ray emission. 

What we have proved to this point is that the bulk Compton 
can account for the X-ray data as long as the jet is pair rich. 
Further constraints about this emission mechanism come from 
the investigation of the pair loading processes occurring in the 
protojet, e.g. photon-photon interaction (Svensson, 1987) or 
the interaction with the X-ray corona field (Sikora & Madejski, 
2002), but this is beyond the goal of this paper 
However a simple viability test of the bulk Compton scenario 
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Figure 8. LECS+MECS+PDS 3 1 -October - 2 November spec- 
trum of BL Lac. fitted by the model in equation|9j using Ti = 1 . 
If Ti < 1 the data are badly fitted by this model. 



can be performed which is based on the X-ray spectral shape. 
If this mechanism is working, in the X-ray band we should 
see the exponential tail of the blueshifted multi-temperature 
blackbody emission of the accretion disk superposed to the 
usual hard inverse Compton emission. Therefore we tried to 
fit BeppoSAX data with a new model described by: 



A(y): 



-Nho-(v) 



-[(vo-v)/vi] 



(9) 



where the second power law parameters are the best fit val- 
ues obtained from modelling PDS data alone with a power law 
model (Fz = 1.76; k2 = 3.19 x 10"^ photons/keV/cm^/s); Nh 
was fixed to 2.5 x 10^' cm"^. This model is not able to repro- 
duce the data unless Fi > 1 (see e.g. in fig. |8]the residuals below 
1 keV). This implies that a blueshifted Shakura-Sunyaev disk 
spectrum cannot reproduce the observed soft X-ray spectrum 
unless being characterized by a soft power law slope in the UV 
band. 



4.3. Two synchrotron components 

We considered then the hypothesis that the optical and the X- 
ray spectra we observed are created by two populations of emit- 
ting electrons possibly located in two regions at different dis- 
tances from the nucleus. One produces a synchrotron compo- 
nent peaking in the IR band and is responsible for the WEBT 
section while the second accounts for the BeppoSAX data. In 
order not to overproduce the optical flux, the X-ray component 
must peak above the UV band. This implies that the emitting 
particle population must have a break (ybreak)- If N(y) oc y^^ 
below 7break, wc require ybieak > L6 x 10^ d^^^B^^^^ to not 
overproduce the optical flux. 

The observed fast X-ray variability implies that the presence 
of a break in the particle distribution could not be an effect 
of the cooling. Since fvar ~ 6 x 10^* s (see table |3 the X-ray 
emitting region must be located at a distance from the nucleus 
R < 6 c fvai/iA ~ 2x10'^ Fjg cm, where ifi ~ 1/F is the jet open- 
ing angle and 6 = 1 /(F( 1 -y6 cos dobs)) is the Doppler factor. The 



dimension of the BLR of BL Lac is Rblr ~ 8.3 x 10'^ Mbh.s 
cm (Bottcher & Bloom, 2000), where Mbh,& is the black hole 
mass in units of 10*^ solar masses. Since Mbh,s ~ 1-7 (Woo & 
Urry, 2002), Rblr ~ 1.4 x 10'^ cm: the X-ray emitting region 
should be inside the Broad Line Region, where the cooling rate 
is high. 

If we assume that a continuous power law distribution of par- 
ticles N(y) oc y " is injected for a time fyai, then an electron 
above jc can cool in a time fyai-. After this time, above the 
particle population steepens to N(y) oc y-'"+i)^ while below y^ 
it remains unchanged. The value of y^. is: 

1 



7 X 10 ^^ 



F?n 6 10 



10 



(10) 



We can 



where we have used f/^^ - 1 x 10"^ Fjg erg cm" 
see that y^ is hardly consistent with the required ybreak, but, if 
during our observation the luminosity of the lines was fainter 
than the average Lblr we have considered, then it is possible 
that -yc ~ ybreak- In this case the required break in the distri- 
bution can be accounted by the radiative cooling. Otherwise, 
not to overproduce the optical flux, the X-ray emitting particle 
distribution should be injected with an intrinsic break. 



4.4. Klein-Nishina effect on tfie synchirotron 
spectrum 

Since the Compton cooling rate of electrons with energies 

m,c2 10^ 



y > yi 



"ThYBEL Fi()VB£ij5 



(11) 



is reduced because of the decline of the Compton scattering 
cross section in the Klein-Nishina regime, for a given electron 
injection rate, the energy distribution of electrons hardens be- 
yond yi . The effect is particularly strong for / = Ug^^^/U'^ » 1 
and causes flattening of the synchrotron spectrum at 



V > Vi 



3.6 X 10^ UzBylS - 



^ 3 ^ ^Qi5 {Ublr/0.005 erg cm-^)'/^(^/r) 



(//lO) 



1/2 v2 



(12) 



The hardened part of the electron energy distribution extends 
up to y2 ~ /'^^yi- At y > y2 the cooling of electrons becomes 
to be dominated by synchrotron radiation and the slope of the 
electron energy distribution becomes the same as at y <s: y\. 
Consequently, at v > V2 ~ fv\ , the synchrotron spectrum steep- 
ens, regaining the slope from the optical band but vertically 
shifted up by a factor / above the extrapolation of the optical 
spectral portion. Such a scenario can reproduce the observed 
spectral "glitch", provided / -10. Details of the model will be 
presented elsewhere; here we would like to mention only that 
a similar scenario has been proposed by Dermer and Atoyan 
(2002) to explain spectral glitches between the optical and X- 
ray spectra of kiloparsec scale jets. 



12 



Ravasio et al.: BeppoSAX and multiwavelength observations of BL Lacertae in 2000 



5. Conclusions 

BL Lac has been the target of a multiwavelength campaign 
during 2000, extending from June to November. BeppoSAX 
observed the source twice, in July and at the end of October, 
while in different optical state of activity. As evidenced also 
from previous observations, BL Lac displays a very complex 
behaviour. We summarize the BeppoSAX results as follows: 

- In July, while optically weak, the source displayed a faint, 
hard X-ray spectrum with positive residuals towards low 
energies. The soft X-ray flux varied in a timescale of a few 
hours, while the hard X-ray flux was almost constant. 

- In October, while BL Lac was bright in the optical band, 
we observed in the X-ray the transition from an extremely 
strong, soft component to a hard component dominating 
above 10-20 keV. The soft spectrum displayed an er- 
ratic temporal behaviour with large and fast variations on 
timescales down to 6 x 10^ s, while the hard component 
remained almost constant. 

- BeppoSAX spectral results are conflrmed by simultaneous 
RXTE data: PCA and MECS spectra are consistent within 
the uncertainties. 

The frequency dependent temporal variability is consistent 
with the spectral analysis results: we observed fast and large 
flux variations when the soft component, interpreted as syn- 
chrotron emission, was dominating, while we observed con- 
stant light curves in the hard section of the spectra which can 
be reproduced as inverse Compton radiation in a leptonic jet 
model. This can be explained since synchrotron X-ray emitting 
electrons are more energetic than those that are producing in- 
verse Compton emission at X-ray energies, so they cool faster. 
This behaviour was already observed in BL Lac (Ravasio et 
al., 2002) as well as in other similar objects, such as ON 231 
(Tagliaferri et al., 2000) or S5 0716+714 (Tagliaferri et al., 
2003). 

The analysis of the multiwavelength Spectral Energy 
Distributions and the comparison with other historical SEDs, 
evidences the exceptionality of the X-ray spectrum of October 
2000: during this observation BL Lac was displaying the high- 
est soft X-ray flux ever recorded and an integrated [2-10] keV 
flux which was only sligthly smaller than that detected in July 
1997, while BL Lac was in an exceptional flaring state and was 
displaying a hard X-ray spectrum (Madejski et al., 1999). 
Moreover, the SED of October 2000 displayed another very 
interesting feature: the soft X-ray data laid above the extrap- 
olation of the optical spectrum, while they should be both 
produced via the same synchrotron emission. To account for 
this inconsistency we have investigated 4 possibiUties, among 
which, however, we cannot at present discriminate: 

1. The dust-to-gas ratio towards BL Lac is higher than the 
interstellar one. However, since this misalignement is not 
seen in the other multiwavelength SED of BL Lac with a 
weU defined synchrotron component (except that of 1980, 
where data have large uncertainties), a sudden large in- 
crease of the dust-to-gas ratio is needed. This hypothesis 
is not unlikely since dusty galactic nebulosities are indeed 
observed towards the source. 



2. This is the first detection of the so called bulk Compton 
emission. In the hypothesis that the jet is pair loaded, we 
can reproduce the observed soft X-ray spectrum via the 
bulk Compton mechanism while the optical and the hard 
X-ray spectra can be explained via the usual synchrotron 
and inverse Compton models. However, a blueshifted 
Shakura-Sunyaev disk spectrum cannot model BeppoSAX 
data. Nevertheless, we cannot exclude this scenario since 
the effective spectral shape of accretion disks in the far UV 
is still uncertain and could have a power law rather than an 
exponential slope (see e.g. Laor et al., 1996). 

3. Two synchrotron components present at different distances 
from the nucleus account for the optical and the X-ray 
spectra, respectively. In order to not overproduce the op- 
tical flux and to account for the X-ray fast variability, the 
X-ray emitting particle population must have a break at en- 
ergies greater than y ~ 1.6 x lO'^S'l'^B'^'^. The break can 
be accounted by the particle cooling or must be intrinsic if 
the line luminosity level is lower or higher than the calcu- 
lated average value, respectively. 

4. Through the synchrotron emitting population, we observe 
the transition from inverse Compton cooling dominated 
to synchrotron cooling dominated particles. Between these 
two conditions, there is a range in which the electrons are 
dominated by an inefficient inverse Compton cooling in 
the Klein-Nishina regime: this produces a hardening in the 
particle population and therefore in the synchrotron spec- 
trum which accounts for the observed optical to X-ray mis- 
aUgnement. 

In the SED of October 2000 of BL Lac, above the synchrotron 
peak, we have detected a well determined misalignement 
between two sections of the synchrotron spectrum. To explain 
this feature, we suggested four scenarios presenting different 
favorable and critical points: at present, however, the theo- 
retical models and the data do not allow us to discriminate 
between them, although the hypothesis of a variable dust-to- 
gas ratio is the most plausible one. 
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